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INPUT 



UNIVERSITY OF ILLINOIS 
DIGITAL COMPUTER 

LIBRARY ROUTINE M 20 - 2*k 
By C. W. Gear 

Eigenvalues of a Symetric Matrix by Givens' Method (SADOI Only) 

Automatic 

507 + Routines N12, Rl, Pl6 and Yl 

128 x 128 

Approximately n /l6o + nf 2.0,000 minutes 

Depends on matrix. So far as known, good to 11 decimal 

places for n less than kO and good to 10 otherwise. 

The matrix must "be scaled so that its norm is less than 1. 

Lj a , . s. j_ 

FF O36 in loc 322 on R.H^S. Sura check failed on read in. 

FF 037 in loc IKK Norm of Matrix has exceeded 1. 

FF O38 in loc 07^ Sum of squares check has failed* (see 
description) 

2k Ol^ on R.H.S. of loc ION. 

Routine is now ready to read in another matrix. 

Normal input the program tape (it uses SADOl). It will 

stop on the word 3k Oil FF O36 in loc 322 (l6), on the 

L.H«S. if the tape has read in correctly, on the R.H.S. 

if there has been a reader error. Black switch the data 

tape in which should be written as follows: 

Consider the matrix in its upper triangular form as 



l ll a 12 a 13 



a 22 a 23 



».«*»«*## 



In 



a 



2n 



* * * . * » 



n-l^n-l n.-J-jja. 



n.n 



Each row should be punched, element by element as a sign 
followed by p decimal digits (p < 12) of the element, and 



at the end of each row the character I should "be punched. 
Thus the data tape will have the following sequence on 
its n signed numbers, N, n-1 signed number, N, «..«,« 
3 signed numbers, N, 2 signed numbers, N, 1 signed number, 
N. The master routine will read in the first row of the 
matrix and plant n, the order of the matrix* It will then 
read in up to each N character, and store this as the next 
row of the triangular part of the matrix. If there are 
more than n-1 numbers in this row it will ignore all those 
after the n-lth, if there are less than n-1, say p, the 
p+lth will be the p+2th of the preceeding row, and so on 
to the n-lth interpreted as the nth of the preceeding row. 
It will stop reading in after it has sensed n N's, 
MODIFICATION This routine uses M12 to input the matrix . If it is desired 

to calculate the matrix elements within the Xlliac prior to 
finding the eigenvalues the routine may be modified in the 
following ways 

Replace N12 by an auxiliary subroutine to calculate the 
elements in the order in which they are read in. N12 is 
entered with the orders 



q 



50 p 
50 q 



26 365F 



It occupies locations 365 through hOJ>. These and the 
additional locations 6 70 through 869 may be used by the 
auxiliary subroutine. At no time must the auxiliary 
routine use other locations (except 0, 1, and 2 ) 
Successive rows of the upper triangular form should be 
calculated and stored in Williams memory locations p 
through 870+n-l. p is 870, 871, 872, ..., 870+n-l^ on 
each successive entry . After the first entry to the 
auxiliary subroutine location 386 must contain 870 + n 
as a left hand address, where n is the order of the 
matrix , 



DESCRIPTION 



If required, library routine Rl is in location 356, Pl6 

in location 277 an( i Yl in location 23- Drum locations 

11720 through 12799 may be used. 

To input the auxiliary subroutine delete the sum check at 

the end of the master routine, i.e. copy up to and including 

the word OOF OOF in location 668, and delete the short last 

section starting with OO8OOK0 Follow this with the auxiliary 

routine with a suitable directive 

Entry should be made to the left hand side of location 17;. 

The read-in is under the control of SADOI° The only floating 

addresses used in the master routine are (Pl6), (N12), (Rl) 

and (Yl)° It will finally stop, as before, on a 2k 01^, 

and if this is obeyed it will re-enter the main routine, 

while the auxiliary routine may specify a new matrix (of 

any order less than 129) • 

For full details see the paper by Wallace Givens on this 

subject in the library files , 

The program falls into two parts, the reduction by a series 

of rotations into Jacobi (bordered diagonal or Tri-diagonal) 

form, and the subsequent calculation of the eigenvalues 

from this form by the use of a binary chopping procedure 

involving the calculation of the signature of the matrix at 

each point of the chop I 

Reduction to Jac obi Formo 

Consider a rotation of the matrix about the (i,j) element 

defined by P' A P. . = A 1 Where P. . = 



c. . 



ij 



>s. . 



IJ 



10 
1 



2 2 
c. . -i-S . . = 1 

IJ IJ 



k 



The only elements affected in this transformation are those 
in the ith and jth rows and columns of A. Since the matrix 
is symetric, and this property is retained under an 
orthogonal rotation, -we need only consider the columns of 
A. 

The elements of the rotated matrix A are thuss 

a = a p, q JL ± \ 

a. = a . = c.oa . + s . . °a . p £ i.j* 
i,p p,i ij pi ij pj * r ,d 

a. = a . = c.oa . + s..»a . p / i»j° 
J,P P,0 ij P<3 ij pi ^ r '° 

a. . = c. .(c. . „a. . + s. . .a. .) + s. .(c. . -a. . + s. ..a. .) 

i,i ij" ij 11 ij ±y ij v ij ij ij jj' 

a . . = a . . = c . . (c . . „a . . + s . . °a . ) + s . . (c . . °a . . - s . . oa _, a ) 
ij Ji ij ij ij ij ii ij ij Jj ij ij' 



Since the trace of a matrix is invariant under an orthogonal 
transformation,, 



a.. = a . . + a . . - a . . 
JJ 11 JJ ii 



The values of c. . and s. . have not yet "been assigned <, In 
ij ij * & 

the Jacobi method, also in the library, they are chosen so 
as to make the (i,j) element zero, However, future rotations 
about other elements alter this, and iteration is necessary . 

Suppose that rotations are made successively about the 

elements (2,3): (2,4); . . . (2,n)s (3»s (3,5) • • -(3,n) : (4,5) . . . 

(n-l,n)s making at each step the elements 

(1,3 )s (l,^)i . . . (l,n)s (2,4)s (2,5); - . . (2,n)s (3,5) - • • (n-2,n) zero 

At any step, elements made zero in previous rotations are 

left zero » " ' 

Thus, ail the elements of the matrix are reduced to zero 

except the diagonal and those adjacent to the diagonal. 



This leaves the matrix 



a ! b l °- 



» 8 -T »Q 



\ a 2 \ 
h 2 a 5 b 5 



0, 



b ,_ a - b 
n»2 n«l n~l 

,0 b _ a 
n«l n 



■with eigenvalues equal to those of the original matrix 

Denote this matrix by Bo 

The eigenvalues of this are the roots of the polynomial 

det(B - ul) = f n (u) = 

From the form of the matrix the polynomial can be defined 
by the sequence 

f (u) = 1 

f x (u) *= a ± ~ u 

f Q (u) = (a c - u)f.(u) - bA n (u) 



1 



f k (u) - (a k - u)f k _ x (u) - ^f (u) 



:d 



.K — tj^i ( sjQi 



If tvo successive f, ! s vanish then f =0. 

i n 

To overcome this, set f. = .- b, / if f. . = f „ = 

k Jt-1 k-1 k«2 

and b, , ^ and f. = a. - u if f. . = b, , = 0, 
K-l ' k k &-± k-1 

Then the number of roots of f (x) =0 greater than u is 



equal to the number of agreements in sign (zero considered 

as positive) in the sequence f n (u), f. (u), f (u)j . . . . .f (u) 

u x d, n 

For accuracy the equations are formed in semi "floating 

point* 

2 
b. is formed and scaled so thats 

2 -pi <=^0 
< K -2 q. < 2 u 

— i u i 

or q = and e. = 79 if b a» 



and a - u la scaled so that; 



»01 X 

a. » u = p.2 » < 8. < 39a o < P-? < 1 



and f (u) = 2" ,ti n 1 



If it_ is set as 0, and v. = jr. - k. the machine 
J i i i=»l 

equations are in the form? 



m, = 2 (p.m. n - 2 Al q. .m. .) if A, < 
i N *i i-1 T.-1 1=2' i 



= v.-8. 
l i 



j./«=iii 



= 2 (2^ p,m .. - q.m. n ) if 4 > 
* -*i 1=1 i 1-2' i — 



Zu = V.+V„ _-€. t 

i l i«l i-1 



where A. = v. ,+5.-€. , 

i l+l i i=l 



S. is found by scaling the quantity in parenthesis to 
lie between 1/2 and 1, and v. is found from the 



appropriate equation 



THE CODE 



The program is read into the memory and stored on the drum 
for future use as temporary storage overwrites much of it. 
The triangular part of the matrix is read in row by row 
and packed on the drum according to the following arrange- 
ment; Consider a track of length 65 words } 



location 


4- 


1 


2560 


a ll 





2625 


a l2 


d.d 



n~6k 



6k 



a 66,66 a 66 3 6l 
a, 



67,67 



66^n 
a 67,n 



2690 



2560 + 
6kx6<3 

2560 + 
65x65 



2560 + 
(n-l)65 



l 13 



L l,65 a 2,6$ ' 



l l,66 a 2,66 



a« 

33 



l,n 2,n 



a 65,6 5 
a 65,66 

a 65,n 



If two locations A and B are such that loc(B) - loc(A) = l(mod 6k) 

then B can be utilised in the minimum access time after A has 

been referred to. 

If p is greater than 65^ then a^ is in location 

2560 + 6kx66 + 65 (p - 65) 3 and a^ is in location 

p - 6k + 2560 c 

Therefore loc(a,v- ) - loc (a^ ) mod 6k 

= - 6kx66 - 65(p«65) + p - 6k (6k) 

= 65 2 - 6k v (6k) 
= (6k + l) 2 (6k) 
= l (6k) 



■= 8 « 

Therefore the rows and coliasas of the matrix can "be 
referred to In minimum access time* 
The sum of the squares of the elements of the matrix 
is calculated as it is read in,, and stored as a check 
against the sum after the reduction to Jacobi form. 

After this reduction, the program will hang up if the 

-39 
difference exceeds 100x2 » 

The eigenvalues are then found from the reduced form, 

by storing the square root R of the sum of squares of the 

original matrix as an upper and lover hound on the absolute 

size of the eigenvalues and adjusting these until they 

converge with information gained about the position of the 

roots by the method outlined. 

Taking P(u) to be the number of eigenvalues greater than 

u, the method can best be seen from a flow chart « 

STORAGE ALLOCATION The program occupies locations 17 through 668 and also on 

the drum. During input of the matrix locations 870 - 997 

are used to store the rows • During reduction to Jacobi 

form locations "jk-2 = 869 and 870 - 997 are used as temporary 

storage for the i+1 and j+lth rows of the matrix which are 

altered in a rotation about the i+1, j+1 element « During 

the calculation of the eigenvalues the a.'s are stored in 

1 

7^2-869, the scaled b. 2 in locations 870 - 997 with their 
scaling factors in 6±k - "jk±. Locations 358 - ^85 hold 
the lower bounds of the eigenvalues, while locations 
k-86 -* 613 hold the upper values . 



DAT E June 27. 1957 RTtl/23/3 9 
PROGRAMED B Y J^P O, G fff- 
APFROVED BYy>?£ ( TtUti ^^ 



lgr 



■? 



Entry 



A 



Al 



Read 

Program 

in. 

Sum 

Check 



"> 



B 



Store 
on drum 



Fail 
FF Stop 



Read 
Program 
from 
drum a 



-* 



<r 



Black 
Swc 



V 



c 



Read first row of 
matrix o Set n 
in program. 



Dl 



FF 
Stop 



<r 



Stop 



N© 



H 



Compute 
Eigens- 
values . 
Print , 
(2) 



Yes 



White switch 
G 



Compare S.S. 
< 100 ,2~ 59 



Gl 



No 



* 



FF stop 



D 



Store on 

drum. 

Sum squares 

< 1? 



—i 



Read JacoM 

form from drum, 

Sum squares . 

2 
Scale "b. 



a 



Yes 



After n 



rows 



D2 



E 



Reduce to 
Jacobi form. 

a) 



Read next 
row 



MAIN FLOW CHART 



10 



i=0 



k = i 



\ 



.3 
set drum 
entries , 



i-i+1 



i=n? 



© 



Set independent 
addresses * 




Yes 



is 



F| FINISH 
Yesj 




11 



a, . - = 

1*1,1+1 1 ? 



No 



a . . , , =0? 

i,j+l 



No 



13 

i .iu. .1 1 urn l"» 



Interchange 
rolss ©i* c 
and a, and 
positions ©f 



X jf -L+.A. 

i * j+i 



19 



Set addr* 
dependent 
on p. 



I 



i 



Yes 
1 10 






I » 



xu 



t = z 2 /2 

v = | /(t + 1/2) t / 
1/2 - 

c * c - (c 2 ~l+s 2 )/2c 

cs-lj t = 



«j=i+X to Ace 



i+1 row in Mem= 



Subrout 

-wggJ 



I §. 



Set j *-de$endeat 



4+1 in Ace 



,3+1 row in Mea< 



Brum 
Subrout 



24 



4+1 in Ace 


^ 


Drum 


3+1 row on 


drum 


Subrout 




1 


25 




j 5! j+1 


"*— © 


i 


p6 




3 -n? 


— © 



Yes 



20 



i '■ , _ = c c a . , 
p, i+1 p,i+l 



+ s 



P>J+1 



■fit . . •* "^ G b 8L_ * 1 "°*a j . t 

Pjj+X ?*J+1 g>**i 



21 



p = p+1 

E 



22 



No 



p a n+l?l Yes 



"•* 



2$. 



Cca 



i+1, i+1 i+1, i+1 



j+l,i+l 



=s a 



i+ljj+1 j+l,j+l 

+e & a ! , , . ' 
l+l, J+l 

*4+l, j+l^i+l, i+1 4 * j+ly j+1 
~® i+1, i+1 



© 



27 



i+1 to Ace = 
— Store- — - 



i+1 row am. drum 



Drum 
Subrout , 



f 



(l) Beductiom to JacoM Form 



1 



GjGl 



11 



=(SoS) 



28 
El/2 



23L, 

u = 



i = 0^P=0 

■-1- 

m Q .- 1-2' 



JO 



-59 



.51 



Set i 

dependent 

addresses 



i = 1+1 
i = n+1? 



-•* 



— _ 21 

Scale a x ~^ 

to p. and 8, 



Yes 



JEas 



^x 



*i-i > ° ? 



*i-i - x 



p, 



=59 



No 



56 



M* - * -1+2 

i-1 



-59 



■**l a J 



= P, 



"} 58 



Y i* 5 I 



-€> 



Yes 



25k. 



m. = ? 

1-2 



-No- 



-22 



l-l 



-No 



la 



* \ = v i-i + V € i-i 



Yes 



kl 



k2 



ml , = 1-2 
i-1 



»59 



,Yes 



i« - 



-i =p i 



^3 



v = 8 
i i 



lf4 



46 



*i-i - 0? 



No 



m! * -1+2 
i 



■59 



■**(?) 



Yes 



.to 



if7 



KL > ? 

1-2 — 



"55" 




m« = 1-2 



»59 



2r„ = p.«i, , - q. , .Bio ^°e i 
1 r i i-1 ^i-l i-2 



Exit it- -» it . 



51 



mi = 

1 



"***($y 



No 
_L_ 



k8> 



A. > 0? 



Yes 



50 



Yes 

52 



r, = 0? 
■» i 



-No- 



-Ai 
2r. = p. cm. *e -q. , .1. ^ 
i ■*! i-1 u i-l 1-2 



Exit jt n -» rt . 
2 



55 




Scale r. 

1/2 < |ml | < 1 







<SK 



A 



v. = £, + 8. 
i i 1 



^1 



(^pV-» V, = Z. + €. , -V. n 

Kj^f^ 1 1 1-1 i°l 




(2 ) Eigenvalues . Part I „ 






56 Yes 



(W> m i-i > 



0? 



No 



57 



-*»} ® i > G? 



Ye? 



58 | 



m, > 0? 



— No 



Yes 



.52 



p * P + I 




61 



6o 



k = j 



Set k 

dependent 
addresses 



74 



J = J + 1 



64 



u - L. > Of 



63 Yes 



■£* 



n 



No. 



Yesg^ 



P-k > 0? 



No 62 



»> k 3 n+1? 



L. = u 

k 

_ ^ 



j68 

1 K^Vc*) 



k t -jfr. — 1 



Yes 



69 



No 67 



u - U k > 0? 



Yes 



Set j 
dependent; 

addresses 



,. 70, 



j = n+l? 



Yes 



Black &w. Stop 



No 







i 71 


Yes 






72 


u. 
J 


- L. -3.2 Jy > 0? 




u = 


(Uj + L.)/2 



















-<D 



No 



ja 



Print 
Eigenvalue, 



(2 ) Eigenvalues . Part II 



~ 13 ~ 
The Brum Subroutine 

This is an independent subroutine which stores the triangular part of a symetric 
matrix on the drum as outlined above « Entry is jaade to read (write) the ith row 
except for the first j elements into (from) Williams memory locations y+j through 
y+n -with the following orders $ 



t 


5(J)0 


y 




50 


«L 


t+x 


26 
00 


k 



75 



store i in 1 
plant q+1 
store y in 2 



No 



W 



with i in the accumulator 
76 

Bead 



77 



Read or write ? 



Store 
"basic 



address of 
read order « 



Write 



19 



80 



Store address of 
"basic write order. 



A> — A» 



y 



85 



y + 65 in "boxes 
93,93 * ,9k and 9k x * 

y+i-1 in 97 and 97 ! 

66 (i-l)+2 in 99 

and 99» * 



— n 




yes 



83 
i > 65?}«*- 



Yes 



84 



1 



66 (i-l) + k in 
95 and 95 x * 



■•» 



65r + i -66 + 2 
in boxes 95 and 95* 

y+r in boxes 9k 
and 9^- r „ 



X = 



81 



♦* 




.3 in loc » 
plant link = 
q + 2 

65i in loc* 3 

65(i-l)+k in boxes 
91 and 91' 

y+k in boxes 90 
and 90' 

y+n in 97,97 ? , 

101 and 101 ! 

y+i-1 in boxes 
93, 93', 9^9^* 
98 and 98 f - 





The Drum Subroutine. Part I, 



- Ik 



15 ( ) 



-*z£L 



I 



£L 



86 11 2560 + ( ) 



I 



92 



increase addresses by 1 



I 



.31 



address 90 > ( ) ? 



No 



Yes 



94 



i 



L5 ( ) 


*■" 


i 


95 


86 11 2560 + ( ) 


i 


96 


increase 9^ ^y 1 an ^ 


95 by 65 


* 


97 


address 9^- > ( ) 



Yes 



9V 



1 



No 



90 1 



92' 



*0 ( } 

ze: 



increase addresses by 1 



SSL 



address 90* > ( ) 



No 



Yes 



1 



351 I 



85 11 2560 + ( ) 



9k 



96 ' 



*o ( ) 



increase 94 * by 1 and 95' by 
65 = 



si 



address 9^ 1 > ( ) 



Yes 



98 



r — V3 
t 



L5 ( 


) 






No 


* 






99 


86 11 2560 + ( 


) 






* 






100 


1 increase address 


es 


"by 


1 


* 






101 


address 98 > ( 


) 






Yes 













No 




99 s 



85 11 2560 + ( ) 



98' 



}^L^L 



100' 



increase addresses by 1 



101 ' 



address 98 » > ( ) 



Yes 



-•*■ 



The Drum Subroutine » Part il 



J LOCATION 


ORDER 
KO 017K 


BOX NO, 


NOTES PAGE 1 




i 


17 


JO 100F B ; Store program on drum* 




50 17F j * 


18 


26 25F | 




02 2688F ! 


19 


00 840F 
26 40^F 


] Exit to Co 

1 j 


20 


50 100F 


I j Head program from drum; \ 




50 20F 






21 


26 23F 
02 2688F 






22 


00 840F 
26 ^-04f 
(11) OOK 




Exit to U<> 


®3 


F5 339F 




Prepare to read JacoM form from tern* 




k2 8k? 


i 1 


6k 


k2 85F . i ■ ' 




k2 72F \ 


65 


Fk 3kW j j 




k2 9W \ j 


66 


te 105F 


■ 


k2 78F : 


67 


^2 92F ! ; 




F5 &%F 


i . 


Do 


k2 I07F 


! 
j 


69 


L5 536F 






ko 71F 






70 


IA 53^F 


j I 




ko 77F 




1 


71 


00 OF 




Read from drum locations 2560, 2625, 2626, 2690, | 




00 OF 




2691, ... | 


72 


22 72F 
kO OF 




! 
| 

i 

i 




| 


i 



M 20 



LOCATION 


1 

ORDER | BOX NO. 


noo!ss mm 2 


73 


L5 71F ' 
F4 3$4F j 


After 130 elements haye "been read reset loop 


74 


40 71F ; 

F4 72F | 

t 


and read 256*2, 2563 » 2628, 2629, • •• to end of 


75 


42 72F | 

LO 547F ! 


matrix . 


7^ 


36.'83F 

26 77F , 


Test for end* i 

.1 

i 


77 


00 OF j 
00 OF 1 


< 


78 


22 ?8F 1 

40 OF 

1 


\ 


79 


L5 TTF j ■ 
f4 334f : 


i 


80 


4o 77F 1 
F5 78F 




i 8l 


42 78F ; 





M 20 



82 

83 

84 

85 
86 
87 
88 

89 



LO 662F 
36 71F 
26 663F 
4l 10F 
4l 11F 

41 12F 
50 OF 
L5 IIP 
74 OF 
L4 10F 
40 10F 
S5 OF 
40 HP 
F5 84F 

42 S4F 
42 8§F 
LO 3^F 



Test for 130 elements* 






Sim sgtmres of a. a s 



LOCATION 



90 
91 
92 

93 
9^ 

95 
96 

91 
98 

99 

100 
101 
102 
103 
104 
105 
106 



ORDER 



32 SkF 
hi 13F 
22 91F 
50 OF 
22 92F 
75 OF 
kO OF 
S5 OF 
kO IF 
IA- 13F 
10 39F 

40 l^F 
L"7 l^F 
hk OF 
lA 12F 
^0 12F 
£5 OF 
*K> 15F 
L3 OF 
36 128F 

41 2F 
50 IF 
L5 OF 
00 IF 
kO OF 
F5 2F 
k2 2F 
L5 OF 
32 101F 
10 IF 

iA 338F 

k) OF 
L5 2F 
FO 335F 



BOX NO 



NOTES 



PAGE 3 



M20 



Square b.'s and scale 



Also sum squares 



LOCATION 



ORDER 



BOX NOo 



NOTSS 



MSB k 



M 20 



10? 
108 

109 
110 
111 
112 

113 

114 
115 
116 

117 

118 

119 
120 
121 
122 
125 



22 


107F 


ko 


OF 


F5 107F 


k2 


107F 


F5 91F 


k2 91F 


k2 


92F 


k2 


105F 


LO 


3^5F 


32 91F 


L5 10F 


Lk 


12F 


iJk 


12F 


LO 351F 


ko 


10F 


L7 


35^-F 


L2 


10F 


32 


ll6F 


FF 


56F 


L5 127F 


kQ 6lkF 


kl 870F 


F5 269F 


k2 


120F 


F5 270F 


k2 


121F 


L5 


353F 


ko 


OF 


LI 


353F 


ko 


OF 


F5 120F 


k2 


120F 


F5 


121F 


k2 


121F 



2 2 

Sum a . + twice sum b, - previotis sxm of 

squares . 



Test to see if it exceeds 100.2 



.39 



Set b = 0, scaling factor = 79 . 



28 



Set upper "bound for eigenvalues = square root 

of sum of 



squares, and lower bound = 
i 


= - square rooto 

i 


1 

! 

] 

i 




i 


... , , m 1 



LOCATION 



ORDER 



BOX NO* 



NO'BS 



PME 5 



M 20 



124 


LO 27^F 




36 120F 


125 


4l 273F 




F5 273F 


126 


ko 272F 




26 131F 


127 


00 OF 




00 79F 


128 


L3 IF 




32 129F 


129 


26 100F 




L5 127F 


130 


kO 2F 




23 105F 


131 


kl 221F 




kl 5F 


132 


4l 222F 




kl 3F 


133 


L5 352F 




kO kF 


134 


F5 222F 




42 222F 


135 


L4 339F 




k2 l42F 


136 


LO 3kXF 




j k2 161F 


137 


42 198F 




; l4 3^if 


138 


1 Lk 34lF 




42 I58F 


139 


| k2 206F 




l LO 3^9F 


140 


j 36 231F 




! L5 273F 



29 



30 



31 



32 



u = 



j - 1 



Const . 



If b. =0 set scaling factor « 79 



i « 



Mp - 1 - 2~ 59 , 



i = 1 



Set i dependent addresses 



i = n+1 ? 



LOCATION 



OEDER 



BOX HO, 



HtffiES 



EACffl 6 



M20 



l4l 
l42 

l44 
145 

i4j 
148 
149 
150 
151 
152 

153 
154 

155 
156 

157 



10 IF 

4o 6f 
43 150F 

L5 OP 
10 IF 
L0 6F 
40 IF 
L3 IF 
36 17TF 
X£ IF 

4o St 
F5 150F 
42 150F 
L7 6f 
00 IF 

36 i46f 
L5 IF 

50 335F 
10 IF 
00 OF 
40 IF 
L5 150F 
LO 275F 
40 2F 
L3 3F 
36 216F 
L3 4f 
32 206F 

51 IF 

74 4f 

10 IF 
40 6F 
S5 OF 
40 11F 



33 



34 
39 



Stale a.-u 



*l-2 - ° ? 



*i~l 



= 0? 



Form q. n »m. „ 
6 and 11 



Store double precision in 



nmmiwwmra i " 



LOCATION 


ORDER 


BOX NO* 


H0TSS PAGE 7 


158 


51 3F 




Form p.. -is. .,. Store double precision in 




74 OF 




7 and 12 


159 


10 w 
4o 7F 






160 


S5 OF 
40 12F 






l6l 


L5 5F 
LO OF 


4o 


A i - v i-i + 8 i + e i-i 


162 


L4 2F 
40 8F 






I63 


43 193F 
L7 8F 






164 


42 226F 

42 227F 




Store as scaling factor 


165 


LO 230F 
36 225F 




A . > 40 ? 

1 =» 


166 


L5 227F 
40 171F 






167 


40 182F 






• 


L5 8F 






168 


36 179F 

L5 229F 


48 


A. > ? j 

1 


169 


46 198F 
LI 12F 




f 


170 


10 39F 
LO 7F 


49 


Calculate r . if A. < 
i 1 


171 


00 OF 
00 OF 






172 


l4 6f 
40 6f 






173 


S5 OF 
L4 11F 




■ ■ 


174 


10 39F 
40 OF 







M 20 



LOCATION 


ORDER 


BOX NO, 


— - • ■ 1 

NOTES PAGE 8 


175 


LI OF 
Lk 6f 






176 


kO 6F 








26 1061, 


■ 


.'■'■■.[ 


ITT 


1$ 2t6f 




'.'■'■" " ■ .'! 


! 


**0 2F 




' 1 


ITS 


J^l IF 




1 


■ 


26 153F 




1 


179 


L5 228F 




.".■'! 




k6 198F 




! 


180 


F5 182F 
42 182F 




1 


181 


50 UP 


50 


Calculate p. if A* > ' 


. 


ia ^p 


• 


] 


<J*ViflC=, 


00 Of 




1 




00 c^ 


: : 


J 
1 


183 


DO lp 




i 




1$ 73T 






3* 


*WP 'IBS? 


1 j 


i 
1 


ji ' Xft5 


1 lo asr 1 




. ■ 

j 


: 


2.6 it^p 






186 


L5 6f 


52 


r ± - ? 


l 


36 223F 






187 


S5 OF 


' 


■ ■ 




^0 of 






188 


L5 6f 
40 IF 


53 


Scale p. 


189 


F5 195F 




■ 




k2 193F 




. 


190 | 


lj if 








00 IF 


, 


. ';> 


191 


kO IF 




• 




56 X89F 







M 20 



LOCATION 
192 

193 
194 

195 
196 

197 
198 

199 
200 
201 
202 
203 
204 
205 
206 
207 
208 



ORDER 



L5 
50 
10 
00 

4o 

50 
L5 

4o 

S5 
40 

L5 
LO 
22 

JLTT 

LO 

4o 

22 

Lk 

ko 

L5 
36 
L5 
36 

F5 
42 
26 

L5 
32 
26 

L3 
32 

L5 
32 

4l 



6f 

OF 
IF 
OF 
3F 
3F 

4f 
3F 
OF 

4f 

193F 

275F 

OF 

OF 

5F 

5F 

201F 

2F 

5F 

3F 

205F 

4f 

134F 

221F 

221F 

134F 

4f 

203F 

134F 

OF 

212F 

3F 

210F 

3F 



BOX NO* 



55 



54 
56 

58 
59 

57 
44 



NOTES 



PAGE 9 



M20 



r ±=\ + € i-i - V: 



m i-i * ° ? 



m i > ° ? 



P a P + 1 



in. > ? 

1 — 



*i=l 



? 



m i„ 2 - ° ? 



LOCATION 


ORDER 


BOX NO. 


NOTES 


PAGE 10 


209 


L5 352F 
kO kF 


^7 


m! - 1-2 " 59 
l 




210 


22 201F 

kl 3F 








211 


LI 352F 
kO kF 


k6 


m£ » -1+2 " 59 




212 


22 201F 










L5 IF 


k2 


*i - P i 




213 


kO kF 








21k 


L5 352F 
kO 3F 


kl 


*i-l " X - 2 " 59 






L5 2F 


k3 


T l" B i 




215 


kQ 5F 
22 201F 






; 


216 


LI kF 
32 217F 


35 


m i.l > ° * 


• 


217 


22 212F 
LI 352F 






i 


218 


40 3F 
LI IF 






' 


219 

220 


^0 ^F 

L5 2F 
40 5F 
22 201F 


56 

38 


v. « 6. 
i i 


• 


221 


00 OF 
00 OF 




P 

• 




222 


00 OF 
00 OF 




Counter for stuna sequence, (i) 




223 

i 


L5 J*F 
kQ 3F 


51 


m i = ° 




22k 


kl kF 










22 201F 








225 


L5 226F 




If scaling factor of A. exceeds kO 


avoid 


1 


22 166F 




subtraction in "boxes ^9 and 50 





M 20 



>■- " " 

I^CAKOH 


OKDEE 


BOX NO. 


IK»S PAGE 11 


' - -*- — " 

226 


*1 335F 








50 OF 






227 


kO OF 
10 OP 




Constants 


228 


00 19®F 








00 OF 






229 


00 200F 

00 or 






230 


00 OF 

00 ^3F 




■ 


231 


L5 272F 


60 


jT* j 




kO 27XF 






232 


L5 27IF 


6l 


Set K dependent addresses. 




lb 269F 






233 


k2 23ff 




■ 




*2 262r 


■ 




23^ 


id 2?3J* 
Ih 2T0F 






235 


k2 263f 








k2 26SW 






! 236 


L5 221F 


63 


? . k>© t 


t 

i 


LO 272^ 






! 237 


32 263F 


<7 


u - \>0 t 




L5 or 






j 238 


FO 273F 






I 


36 262F 




! _ ' . ■ ■: 


! 239 


10 IF I 


; "' ..-..■ 


[ 
j 


F5 271F j 6a 


& m k + X 


2*t0 


kO 273J* 




■ 


j 


L0 3W 


62 


k»a + 1 ? 


| 2kL 


36 232F 






1 
1 


15 27^f 






2**2 


iA 27OF 


69 


Set j de$»end$»t addresses 




k2 2^6f 


: 





M20 



LOCATION 



ORDER I BOX NO. 



2k3 


42 257F 




L5 272F 


2kk 


IX 269F 




k2 2^5F 


2k$ 


k2 259F 




L5 OF 


2k6 


LO 355F 




LO OF 


2^7 


32 257F 




92 131F 


2kS 


92 515F 




92 967F 


2^9 


L5 272F 




22 250F 


250 


J2 3F 




50 250F 


251 


26 27T# 




92 6JF 


252 


92 963F 




L5 275F 


253 


50 12? 




50 25JF 


254 


26 277F 




F5 2°F2F 


255 


kO 272F 




LO 5^8F 


256 


52 2ta 




22 268F 


257 


00 OF 




L5 OF 


258 


10 IF 




kO OF 


259 


22 259F 




L5 OF 



71 



73 



N0EB5 



HUB12 



M20 



7^ 
.70 



72 



U. - L. - 3»2~ 59 > t 



Ertat eigenvalues aM. mmfoer of root (j)» 



j = 4 + 1 
J = n + 1 ? 



u = (U. + L.)/2 



LOCATXQI 



260 
26l 
262 
263 

2.6k 
265 

266 

267 
268 
269 
270 
271 
272 
273 
27^ 
275 
276 



10 LP 
L4 OF 
kO 273F 
26 131F 
L5 273F 
40 OF 
22 239^ 
U OF 
IA 273F 
36 266F 
22 239F 

00 of 

X5 273F 
k) OF 
22 239F 
00 OF 
22 268F 
2k 20F 
00 OF 
00 48SF 
00 OF 

00 ym 
00 OF 
00 OF 
00 OF 
00 OF 
00 OF 

00 m 

Ttf 353F 
40 486f 
10 IF 
00 2F 
00 OF 
00 79F 



BOX HO. 



66 



& 



65 



mm 13 



M 20 



u 



u=L>Ot 



m U 



Useful constants 



mmmmmmw mammt^t 



Locmtsi 



333 
33^ 
335 
33^ 
337 



539 
3*K> 

3*U 

3^ 
3^3 

3^3 
3^ 
3^7 
3if8 



(106) QOK 

00 OP 
00 OF 
00 Of 
00 &$F 
0© Of 

oo or 

65 X1F 



L 



00 2560F 

ao of 

00 OF 

00 OF 

00 7to 
30 ffi 

00 OF 

00 IfiflF 

00 128F 
00 OF 

00 OF 
80 OF 
0© OF 

%0 2F 

30 €F 

IS 92^ 

50 OF 

m J2F 

50 OF 
22 72F 
i*0 OF 
BO OF 
00 OF 



i « M)i» ij ii<i M ! OJ# * w a w ^Wj p i i ' m»i i Mwa BW ffl BltBawft* 
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M20 



More useful constants 



ilWiipi^'WWI!!^^ 



i ! .! ^iiurjnirjffi»i»pnmrtii|iwMii W W 



L0GA.TI0H 


ORDM 


BOX 10* 


norm mm 15 


3^9 


00 128F 




'■■■'■ 




00 OF 




■. ■ ':'.'.,■'■■' 


350 


50 OF 
26 k2kl 






351 


00 OF 








00 OF 




' ' : ■ ' "■" ' ' " ' "' ■ ■ ■" -. ■'. ' ' 


352 


7L 4095F 
LL k)95F 




1*2 ~ 39 


353 


00 OF 
00 OF 






35^ 


00 OF 
00 10GF 






355 


00 OF 

00 w 






356 


(11) OOK 






365 


(112) 001 






k& 


kl 9F 

kl 10F 




^ogran starts here 


k05 


kl kf 








22 kC&F 






ko6 


50 8?0F 
50 koSw 






lK)T 


26 365F 




Scad ia the fint row lata ltfeatiewi 




L5 366? 




S70 - 870-m^l* 


km 


k6 350F 

LO 3ta 




S70 + ais location (112) + 21* 


k09 


k6 3m 




Set the test coastaat* depemd«Kt ea a« 




10 20F 




■ 


klQ 


k2 jk6r 






kn 


1* 3ta 


. 






te 3^9F 







M20 



LOCATION 


ORDER 


BOX HO, 


NOTES PAGE 16 


4l2 


42 545F 








LO 34lF 




• 


4l3 


LO 339F 
42 548f 






4l4 


FO 335F 

42 343F 






415 


.42 333F 
FO 335F 






kl6 


42 542F 








L5 446f 




aad. arrange for prograM to skip this is future; 


klf 


4o 407F 








15 4F 




. 


4l8 


42 420F 
1*5 4F 






4l9 


jo 8701 




Store us lag drum siteoutime 




50 419F 






420 


26 555F 
00 OF 






421 


L5 447F 

4o 424f 




- 


422 


L5 4o6f 

46 423F 






423 


50 OF 




Stua squares r takisg tsrlee tbs stas ©f t&e off 




k6 h2k¥ 




diagonal eleaests 


424 


75 OF 

00 OF 






425 


I* 9F 
40 9F 






426 


36 427F 




Stop If sons exceectol. 




FF 55P 






427 


S5 OF 








L-4 10F 






428 


10 39P 








40 GF 







M20 



LOCATION 


OKBEB 


BOX HO, 


boob mm it 


429 


LT OF 








L4 9F 




' . • 


430 


40 9F 








S3 OF 




■ ; -■'■' ■ ■ . 


4 3 1 


40 IGF 




. ■ 




L5 W&F 






432 


40 424f 








1$ 423F 




■ 


433 


L4 44^f 




: 




4^423F 


■ 


' ' ' ■ ■ '' 


434 


W> 35GF 




' ' ' ' . . ■' ■ .'"' ' ■ ■ 




36 43^F 




'■'.- 


435 


15 42^ 




■ ■.-■ ''.■' 




B€ k23$ 


; 


'''■'. : '■■. ' 


436 


40 10F 


i 


' ' . . ' .'' ■ ... '■ .'■■ ■". '; 




m ^6f 




Xaermsre mte«*§** *** tett« 


437 


IA 44^F 




. . '■''.'..■ ' : 




46 4^F 




■ ■ i 


43$ 


F3 4F 




' ■ .. ' .■ ■. ; ... , ' . 




te 4? 


; 




439 


LO 343F 


• 






32 4o6F 




. . ■ ■ . . 


44o 


L5 443F 

4o 407F 




Rese'fc r$&t£jM5 


44l 


L5 444r 


. 




.. 


I0^f 




■■' ■■ ...■■.'••■ '■ 


kk2 


kim 








26 44^ 


■ 


. ' '.■'■' .... '. ■ 

i 


443 


26 36^ 




C@K@tatits 




L5 3B6F 






444 


50 8T«r 








50 4©6F 

■ 




. ' , ; . ■ . . '■ 


445 


a© if 




■ ■ .' ■■ ' . .. '.' . ■ ■. ■■ '' " ■' , ■.' 




©0 ©$• 




■,' .'.■■■ 

. ■■■■•■■ 
• 



M 20 



LQC&TION 


QRDEft 


1CK W)» 


hots mm xS 


446 


26 365F 








22 417F 




■ ■ . 


447 


75 of 

26 425F 






448 


75 OF 
00 IF 






44o 

"fry 


W 9F 
50 44$r 


■ 


Compute square r«©t of siw or squares. 


450 


26 35^F 
40 353F 




Jfcegfei j^mctiom 


451 

4^2 
453 


L5 9F 
40 351F 
43 554F 
L5 554F 
42 468? 


1 


1*0 

k « i 


3 


Set di*ipt si$&r©>ttiB@ salaries 


454 

^55 


k2 475F 
42 543F 
42 5^BF 

15 554F 

42 554F 


4 


i m i+1 


456 


LO 342F 
36 63F 


5 


i « a f 

. ' ■ ■ 


457 
45S 


L5 339F 
1*4 55^F 
42 476F 
L4 341F 
42 477F 
42 344F 


6 


Set i dependent addresses 

: . 


460 


00 20F 






461 


46 54lF 

F5 476F 
00 20F 






462 

> 


k6 532F 

46 535F 

■ 




i 



M20 



tockTicm 


owmr 


I0£ 10^ 




k63 


46 540F 








" 46 516F 




■ ■ ' ■ ''■■■ .' ■ ,'■ . ■ 


k& 


L4 541F 
46 S39f 




■ . . ■ '"■■.' .■;'.; 


k65 


l»6 556F 








46 513T 






k66 


^ 554f 


7 


j » 1+1 




4g 55^r 




. ■ ■ ■ ''■'.''. 


k6j 


50 ^BF 




■ 




50 46?F 




' 


W8 


26 5OT 
©0 OP 




Inter tam jnfrriKtljte t# rmM. ■ i+JL tit 'rwr.- 


*6q 


15 55ST 
L4 5^F 


a 


Set j &e$«»&«ftt address 


Vfo 


42 5l4f 




■ 




te 53W 




: ■ . ■ .' ■ 


*-7X 


4g 536r 




■ ■ ■ ■ . ■ . . 




L4 34U 


■ ■ 


■ .'■'■" 


47s 


te 515F 




■ ■ 




42 5371' 




■ ■■■■■ 


*73 


42 9Mar 








F5 55ftT 


9 


SWier toai s*3&r#iffclae t# r®a& (j+l-)tli t&w 


\1K 


50 $7C3r 








50 V|fcr 




; 


VRS 


26 55^ 




■ ■ ■ • 




00 OOP 




'' ■ ■ ' ' ' ' ' '''■■■.■'' 


*76 


4i of 








ia or 






%7T 


40 -jr 

15 OF 






V* 


4o4f 




■■■ -'' ' '■■■■'■■' : ■ 




13 4f 


ii 


v* 


*W 


36 544F 






LT 3F 


10 


l% i+ xl 2 Krfrtl * 



H£0 



; location 


oisons 


18K SO 


WVTOPC! VA/Hf Oft 


480 


12 4f 






\ 

! . 


36 506F 






481 


50 3? 


13 


Interchange roles of c «kL amd positions 




L5 4F 




of a i, i+ i "* a i,j + i 


! 482 


4o 3F 




■ . ' 


s 


S5 OF 




■ ■■■,■■ 


I 483 


4o 4f 






| 


L5 553J 1 




• 


! km 


42 497F 




■ ■• ' ■ 


1 


46 504f 




i . ' ' ■ ■ 


J 48$ 


50 335F 


14 


C*leul*ie c aad • ■....' 


i 


15 4f 






s 486 


66 ^r 






i 


S5 OP 


■ 


■ 


; m 


40 5? 


■ 






TJ !5F 


. 


' ■ '■ ■ '" ' 


488 


10 XT 




■■'.'■■ ■ ■• . 




40 IF 




■ ' .' ' 


48$ 


L3 IF 






I 


36 505F 






! 490 


U IF 






i 


40 IF 




■ 


491 


50 338® 




' ' . 


1 


6h W 




-■■■. ' • ' ' '■ 


49S 


m of 

36 494F 






493 


26 5<^ 




■ ■ 




00 Of 






494 


S5 Of 






! 


50 494F 






^ 


26 356F 








40 IF 




. . 


>** 


50 IF 






| 


7J 5F . 






| 







M20 



LOCATION 



497 
498 
499 
500 
501 
502 
503 

505 



507 
508 

509 
510 

513 



ORDER 



OF 
OP 
OF 
OF 
OF 
IF 



BOX HO, 



ko 
ko 

50 

79 
40 
50 
79 
LO 

I* 
50 
66 

S5 
10 

ih 
ko 

22 

&5 
22 
kQ 

L3 
56 
26 

L5 
26 

L5 



4-0 OF 
L5 553F j 
ho 5F j 
40 6f ! 
L3 OF 
32 513F 
22 514F j 
LI 553FJ 



OF 
335F1 
IF 
OF 

IF 

IF 

OF 

51^F 

33% 

495F 

IF 

IF 

509F 

508F 

550F 



15 



NOTES 



PA0E 21 



M 20 



16 



17 



■*"jf*J**" ■ 



i,i+l 



I ? 



l i,j+l " a i,i+l ? 



C a s = 2 2 



c «'*■ a 2 2 



LOCATION 
51k 

515 
516 

517 
518 

519 
520 

521 

522 

527 
528 
529 
530 



ORDM 

40 6F 

L5 OF 
40 OF 

13 or 

1* OF 
40 Of 
L5 3Wtf 

42 526f 
42 530T 

42 524F 
LO 54lF 

42 ^sr- 
0020F 
46 528F 
46 523F 
LO 540F 
56 532F 
26 523F 
50 OF 
TJ 6F 
40 IF 

50 or 
7JT9T 

LD.1F 

40 IF 
50 0? 

40 2F 
50 OF 
7J5F 
L4 2F 
40 OF 

L5 IF 
40 OF 



BOX NQ< 



18 
If 



22 



20 



HOTES 



PAGE 22 



M20 



p - i 
Set p 



addresses 



p = m+1 % 



l©t&te mtvis, s&atiplisrbim l©a$ 





ORDER 


bck io. 


IOTES 


531 


f5 5S*€F 
22 517? 


21 


*«#*+! 


532 


50 OF 


23 






TJ 5F 




1 

1 


535 


4o u 




■ ! 

■ ! 




50 Of 






53^ 


l4 if 






535 


4o or 




' 




50 5P 


■ '■ 




536 


7J0F 
40 IF 






53? 


50 6F 








TJ OF 






338 


lA IF 








40 OF 




■■■■'''. . '■ 


539 


40 OF 
15 OF 


' 


■ i 

■"■■■■■ ! 


540 


w m 

40 OF 




■ ■ ■■'■■'! 


5*1 


41 OF 








F5 552F 


24 


Store 4+1 %&. row m drm 


542 1 JO B7GF 

1 






1 50 5*»-2F 






5*3 


26 555F 

00 OF 




■ , ■ ' ' ' : 


5* 


F5 5W 


25 


5 « >1 




42 532F 


. 




5*5 


LO 543F 


26 


j»B! 




36 4€9F 






5W 


F5 554F 
22 54TF 




■ ■ . 


5+7 


JO 7**2F 


27 


itore^iH-J^tftuQiL . dtfiai 




50 54?F 







MGE 23 M 2C 
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LOCATION 
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5^8 


26 555F 








00 OF 






5^9 


22 ^52F 
00 OF 






550 


00 5F 




Ccmsta&ts 




00 6f 






551 


00 6f 
00 5F 






552 


00 OF 

00 OF 




. 


553 


5K 2087F 








99 2557F 






55^ 


00 OF 








00 OF 




' ! 


555 


kO IF 




Start of drum subroutine 




K5 5F 




■ 


556 


k2 56^F 
10 20F 


75 


Store i, q+1 and y 


557 

■ 


42 2F 
32 658F 


76 


Read or write ? 


558 


F5 624p 


77 


Store address of basic read order 




k2 568F 






559 


k2 589F 




■ 




4-2 596F 




\' 


560 


42 607F 








L5 623F 


78 


K = X ? 


561 


42 618F 








42 6lOF 






562 


L4 555F 








42 621F 




' 


563 


F4 555F 








42 £20F 






564 


4l OF 


81 


Form addresses for read (-write) loops 




L5 OF ; 
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ORDER 
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565 


k2 OF 
F5 564f 






566 


k2 64lF 
50 IF 






567 


75 33^F 
S5 636F 






568 


ij-0 3F 

L5 OF 






569 


LO 33^F 








lA OF 






570 


lA 3F 








40 625F 






571 


40 643F 








iA 2F 






572 


iA OF 
42 626F 






573 


k2 642F 






: 


L5 533F 






57^ 


iA 2F 








42 613F 






575 


42 6l4F 
42 615F 






576 


k2 616F 




. ' . 




L5 2F 




, ' ■ 


577 


iA IF 








FO 335F 




■ 


578 


k2 6llF 
42 612F 




■ 


579 


42 63IF 
42 61+7F 






580 


42 637F 








42 653F 







M 20 



LOCATION 



ORDER 



BOX NO. 



NOTES 



PA$S &£ 



M20 



581 


Ll IP 




L4 334F 


582 


32 605F 




L5 2F 


583 


L4 334f 




42 6llF 


584 


42 6l2F 




42 63XF 


585 


42 64?F 




L5 637F 


586 


42 6lj5F 




42 6l4F 


587 


1*5 620F 




42 65OF 


588 


L5 567F 




42 633F 


589 


L5 3F 




L4 OF 


590 


L4 IF 




LO 6X9F 


591 


40 636F 




4o 6^4f 


59S ■ 


L5 OF 




LO 334F 


593 


56 603F 




L5 618F 


594 


42 602F 




L5 617F 


595 


46 597F 




46 601F 


596 


50 334F 




L5 OF 


597 


74 OF 




£15 64lF 



83 



85 



86 



88 



S9 



i > 65 ? 



Form addresses for reeld (wlte) loops 



k> 65 



T * 65 



Construct addresses 
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ORDER 
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M 20 



598 

599 

600 

601 

602 

603 

60V 

605 

606 

607 
608 
609 
610 
611 
612 
613 
61k 



Lk 338F 
IA IF 
LO 667F 
kO 63OF 
kO 6k8F 
L5 2F 
lA OF 
k2 63IF 
k2 6kjF 
26 OF 
L5 621F 
k2 602F 
M 597F 
^7 601F 
26 596F 

L5 597F 
^2 633F 
k2 65OF 
L5 3F 
iA of 
ih if 
FO 33^F 
ko 630F 
*K) 6^8f 
22 6lOF 
26 OF 
K2 626F 
kO OF 
K2 6^-2F 
L5 OF 
22 63IF 
kO OF 
22 6^7F 
L5 OF 



87 



■■>!).* r " ^ 



*j "-^ 



81+ 



4 S ^ 



Test constants used in boxes 93, 97, 101, 
93*, 97', and 101 T . 



LOCATIOH 


ORDER 
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615 


K2 637P 
ko OP 






616 


K2 653F 








L5 OF 




■ 


617 


00 35^F 




Constant* 




00 OF 






618 


00 OP 








00 OF 




■ ■ 


619 


00 OF 








01 258F 






620 


00 OF 








00 OF 




'.■■'.' 


621 


00 OF 




■ '.■''■ 




00 OF 




■ ' ■ • • ; ■ ■ ■ ■ ' 


622 


00 OF 

00 625F 






623 


00 OF 
00 6k2F 






62k 


00 OF 




■ ''.■'■.' 




00 336F 




. ■ . ■ ••■ ^ . ■ 


625 


00 OF 
00 OF 


91* 




626 


22 626F 


90' 






kO OF 






627 


F5 625F 


92' 






ko 625F 




■ ■ 


628 


F5 626F 




Head loops 




k2 626F 






629 


LO 611F 


93' 


* 




36 625F 






630 


00 OF 

00 OF 


95' 


■ 


631 


22 63IF 


9k ! 






kO OF 







M 20 
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632 


F5 63IF 


96 « 


■ . . 


■ 


42 631F 




.■■'■'" l 


633- 


LO 613F 


97' 






36 OF 


. 




63k 


L5 63OF 




:■ ' '' ■ . ' - . ..' '' 




L4 334f 






635 


40 630F 




■ .. • .-.'' 




26 63OF 




.. ■ ■ . . 


636 


00 OF 


99 1 






00 OF 






637 


22 637F 


98' 


. ■ . ■ 




40 OF 


■■ 




638 


F5 636F 


10©' 


■ "■' - '■ .. ■'■ - ." . 


, 


40 636F 


' ■ . ■ 




639 


F3 $37F 


.' 


■ ' ■ ■ . 




42 63TF 


! 




6*1-0 


LO 6l5F 


101. , 


; " ' ' ' "' : 


1 


36 636F 


; 


' 


64i 


22 641F 


102 






26 OF 






642 


22 642F 


90 


: 




L5 OF 






643 


00 OF 


91 






00 OF 




: ' 


644 


F5 643F 


92 


Write leaps 


' 


40 643F 




'■■'■'•■ .; : : 


645 


F5 642F 




: ■ . 

1 




42 642F 




1 ■ ■ . . ' ; 


646 


LO 612F , 


93 


. ' ■ ■ ■■■ . 

• 




32 642F : 




! ' . ' 


647 


22 64?F 
L5 OF 


94 


■' '■■■'■ 


648 


00 OF 


95 






00 OF 







M 20 



LOCATION 



BOX HQ« 



NOUS 

i iiri i in i i i i.i i WWi iii i iwii i i im ii i n i mmm mmnMW 



fJJM 50 H 20 

-- ■-•-»■-■ — -p- i i. i -hi i . 



6k9 
650 
651 
652 
653 

655 

656 

657 

658 

659 
660 
661 
662 
663 
66k 
66^ 



F5 64-7F 
42 6kjF 
LO 6l4F 
32 OF 
L5 648F 
L4 33^F 
4o 6kBF 
22 647F 
22 6^3F 
V? OF 
00 OF 
00 OF 
F5 6^4f 
40 65%F 
F5 653F 
42 653F 
LO 6l6F 
32 653F 
22 €4'1F 
L5 6s4f 
42 568F 
42 589F 
42 596F 
42 60?F 
L5 622F 
26 56IF 
IS 78F 

4o 936F 

F5 336F 

F4 3J55F 

40 71F 

F5 71F 
40 77F 
L5 668F 



5* 

97 



98 

99 
100 



101 



79 



80 



Store address of "basic vrite order 



Test constant ueed in read "back of Jacobi form 

Reset same after 130 elements have "been read, 
and adjust test const so this block is disobeyed 
In future. 



LOCATION 


ORDER 


BOX NOc 


HOTES 


PAQS 51 


666 


kO 662F 

26 ?1F 








667 


00 F 

01 193F 








668 


00 F 

00 F 










00 800& 




Sum check on Routine . If It Is 
modified, delete from here on. 


to "be 


800 


81 40F 
12' T 








801 


kO F 

L3 F 




; 


"" ■ 


802 


3k 17F 

FF 5^F 

26 800N 








■ 


58 ?5^3NL 


^N 


Sum check const. 
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